In [3]:
from pygadgetreader import *
import numpy as np
import pandas as pd
#import matplotlib
#matplotlib.use('Qt5Agg')
import matplotlib.pyplot as plt
import time
%matplotlib inline
import healpy as hp
from sklearn.neighbors import KDTree
In [4]:
# all scales are in Mpc/h
stackRadius = 2
min_box = 0.0
max_box = 500.0
NSIDE = 32
min_edge = stackRadius*np.ones(3)
max_edge = max_box-stackRadius*np.ones(3)
1. Read halos¶
- and extract coordinates. Here we may restrict it to certain halo size/mass range
In [5]:
dataDir = '/mnt/net/novo2/data/astro/nbody/millennium/'
haloFile = dataDir + 'halos.csv'
halos = pd.read_csv(haloFile,sep=',')
h_xyz0 = halos[['x','y','z']].values
2. Read Millenium and build tree¶
- read the given block
- determine bounding box
- make kd-tree
- collect neighbors of halos
- calculate
In [6]:
def healPixStat(blockId):
infile = '/mnt/local/data/astro/nbody/millennium/snapdir_063/snap_millennium_063.' + str(blockId)
coordinates = readsnap(infile, 'pos', 'dm', single=1) ## single=1 makes it read only one block
min_xyz = coordinates.min(axis=0)
max_xyz = coordinates.max(axis=0)
# generate kd-tree
tree = KDTree(coordinates, leaf_size=2)
# select only halos in this block
inidx = np.all((h_xyz0 >= min_xyz-stackRadius) & (h_xyz0 <= max_xyz+stackRadius) & (h_xyz0>min_edge) & (h_xyz0<max_edge) , axis=1)
h_xyz = halos[inidx][['x','y','z']].values
# get points around halos
ind = tree.query_radius(h_xyz, r=stackRadius)
# relative vectors
neighvec= []
for i in range(ind.shape[0]):
xyz = h_xyz[i]
neighvec.append(coordinates[ind[i]]-xyz)
neighvec = np.vstack(np.array(neighvec))
# healpix binned counts
healIdx = hp.vec2pix(NSIDE, neighvec[:,0],neighvec[:,1],neighvec[:,2])
healIdxList, healCounts = np.unique(healIdx, return_counts=True)
return healCounts
In [155]:
def healPixStatHalfMassR(blockId):
infile = '/mnt/local/data/astro/nbody/millennium/snapdir_063/snap_millennium_063.' + str(blockId)
coordinates = readsnap(infile, 'pos', 'dm', single=1) ## single=1 makes it read only one block
min_xyz = coordinates.min(axis=0)
max_xyz = coordinates.max(axis=0)
# generate kd-tree
tree = KDTree(coordinates, leaf_size=2)
# select only halos in this block
inidx = np.all((h_xyz0 >= min_xyz-stackRadius) & (h_xyz0 <= max_xyz+stackRadius) & (h_xyz0>min_edge) & (h_xyz0<max_edge) , axis=1)
h_xyz = halos[inidx][['x','y','z']].values
h_rad = 4*halos[inidx]['halfMassRadius'].values
# get points around halos
ind = tree.query_radius(h_xyz, r=h_rad)
# relative vectors
neighvec= []
for i in range(ind.shape[0]):
xyz = h_xyz[i]
neighvec.append(coordinates[ind[i]]-xyz)
neighvec = np.vstack(np.array(neighvec))
# healpix binned counts
healIdx = hp.vec2pix(NSIDE, neighvec[:,0],neighvec[:,1],neighvec[:,2])
healIdxList, healCounts = np.unique(healIdx, return_counts=True)
return healCounts
In [118]:
#for i in np.arange(30,dtype=int):
#for i in np.arange(30,dtype=int)+30:
#for i in np.arange(30,dtype=int)+60:
#for i in np.arange(30,dtype=int)+90:
#for i in np.arange(8,dtype=int)+120:
#for i in np.arange(32,dtype=int)+128:
for i in np.arange(352,dtype=int)+160:
print(i,",",end='')
start = time.time()
if i == 0:
stack = healPixStat(i)
else:
stack += healPixStat(i)
end = time.time()
print("%f s." %(end-start))
160 ,Returning DM Positions 40.376148 s. 161 ,Returning DM Positions 61.177875 s. 162 ,Returning DM Positions 47.030558 s. 163 ,Returning DM Positions 48.714451 s. 164 ,Returning DM Positions 46.445572 s. 165 ,Returning DM Positions 46.154755 s. 166 ,Returning DM Positions 45.393591 s. 167 ,Returning DM Positions 47.972569 s. 168 ,Returning DM Positions 40.132210 s. 169 ,Returning DM Positions 39.327806 s. 170 ,Returning DM Positions 51.713339 s. 171 ,Returning DM Positions 51.616596 s. 172 ,Returning DM Positions 46.106400 s. 173 ,Returning DM Positions 43.549723 s. 174 ,Returning DM Positions 36.891017 s. 175 ,Returning DM Positions 164.785084 s. 176 ,Returning DM Positions 40.639476 s. 177 ,Returning DM Positions 39.988624 s. 178 ,Returning DM Positions 39.720730 s. 179 ,Returning DM Positions 57.232196 s. 180 ,Returning DM Positions 46.783262 s. 181 ,Returning DM Positions 57.582483 s. 182 ,Returning DM Positions 47.477057 s. 183 ,Returning DM Positions 38.232492 s. 184 ,Returning DM Positions 48.927846 s. 185 ,Returning DM Positions 50.019899 s. 186 ,Returning DM Positions 49.801412 s. 187 ,Returning DM Positions 40.833246 s. 188 ,Returning DM Positions 47.392954 s. 189 ,Returning DM Positions 45.745402 s. 190 ,Returning DM Positions 59.497066 s. 191 ,Returning DM Positions 40.516515 s. 192 ,Returning DM Positions 44.267869 s. 193 ,Returning DM Positions 45.102104 s. 194 ,Returning DM Positions 47.297786 s. 195 ,Returning DM Positions 40.949105 s. 196 ,Returning DM Positions 56.835732 s. 197 ,Returning DM Positions 60.764197 s. 198 ,Returning DM Positions 46.969229 s. 199 ,Returning DM Positions 229.545874 s. 200 ,Returning DM Positions 48.124481 s. 201 ,Returning DM Positions 58.046249 s. 202 ,Returning DM Positions 48.898693 s. 203 ,Returning DM Positions 45.534227 s. 204 ,Returning DM Positions 54.266915 s. 205 ,Returning DM Positions 60.980572 s. 206 ,Returning DM Positions 45.646584 s. 207 ,Returning DM Positions 38.600029 s. 208 ,Returning DM Positions 46.229537 s. 209 ,Returning DM Positions 57.627896 s. 210 ,Returning DM Positions 168.011387 s. 211 ,Returning DM Positions 67.264925 s. 212 ,Returning DM Positions 48.788373 s. 213 ,Returning DM Positions 54.088480 s. 214 ,Returning DM Positions 39.447329 s. 215 ,Returning DM Positions 60.775678 s. 216 ,Returning DM Positions 61.827881 s. 217 ,Returning DM Positions 65.170526 s. 218 ,Returning DM Positions 54.167471 s. 219 ,Returning DM Positions 45.048825 s. 220 ,Returning DM Positions 52.320724 s. 221 ,Returning DM Positions 46.964803 s. 222 ,Returning DM Positions 330.223938 s. 223 ,Returning DM Positions 44.015862 s. 224 ,Returning DM Positions 54.671062 s. 225 ,Returning DM Positions 47.199643 s. 226 ,Returning DM Positions 54.234487 s. 227 ,Returning DM Positions 44.855694 s. 228 ,Returning DM Positions 45.339900 s. 229 ,Returning DM Positions 52.030188 s. 230 ,Returning DM Positions 64.381773 s. 231 ,Returning DM Positions 594.140021 s. 232 ,Returning DM Positions 45.926254 s. 233 ,Returning DM Positions 43.143049 s. 234 ,Returning DM Positions 46.629971 s. 235 ,Returning DM Positions 38.580859 s. 236 ,Returning DM Positions 48.132576 s. 237 ,Returning DM Positions 42.068217 s. 238 ,Returning DM Positions 58.407591 s. 239 ,Returning DM Positions 65.970589 s. 240 ,Returning DM Positions 45.127771 s. 241 ,Returning DM Positions 58.430055 s. 242 ,Returning DM Positions 59.325354 s. 243 ,Returning DM Positions 113.792433 s. 244 ,Returning DM Positions 55.177770 s. 245 ,Returning DM Positions 43.958232 s. 246 ,Returning DM Positions 86.502871 s. 247 ,Returning DM Positions 61.440560 s. 248 ,Returning DM Positions 47.253617 s. 249 ,Returning DM Positions 132.216079 s. 250 ,Returning DM Positions 52.092284 s. 251 ,Returning DM Positions 63.257602 s. 252 ,Returning DM Positions 44.734559 s. 253 ,Returning DM Positions 50.543836 s. 254 ,Returning DM Positions 42.283497 s. 255 ,Returning DM Positions 267.392662 s. 256 ,Returning DM Positions 48.263378 s. 257 ,Returning DM Positions 52.342079 s. 258 ,Returning DM Positions 44.068656 s. 259 ,Returning DM Positions 73.248613 s. 260 ,Returning DM Positions 48.576206 s. 261 ,Returning DM Positions 39.254416 s. 262 ,Returning DM Positions 41.551976 s. 263 ,Returning DM Positions 48.536269 s. 264 ,Returning DM Positions 328.963822 s. 265 ,Returning DM Positions 56.950091 s. 266 ,Returning DM Positions 42.372515 s. 267 ,Returning DM Positions 53.429865 s. 268 ,Returning DM Positions 45.078347 s. 269 ,Returning DM Positions 47.988992 s. 270 ,Returning DM Positions 63.957458 s. 271 ,Returning DM Positions 51.271504 s. 272 ,Returning DM Positions 45.340841 s. 273 ,Returning DM Positions 57.978406 s. 274 ,Returning DM Positions 42.334137 s. 275 ,Returning DM Positions 48.797336 s. 276 ,Returning DM Positions 45.132507 s. 277 ,Returning DM Positions 51.408332 s. 278 ,Returning DM Positions 66.696973 s. 279 ,Returning DM Positions 341.478273 s. 280 ,Returning DM Positions 43.742980 s. 281 ,Returning DM Positions 50.621770 s. 282 ,Returning DM Positions 42.363271 s. 283 ,Returning DM Positions 93.108260 s. 284 ,Returning DM Positions 46.436381 s. 285 ,Returning DM Positions 45.607444 s. 286 ,Returning DM Positions 46.278139 s. 287 ,Returning DM Positions 54.612319 s. 288 ,Returning DM Positions 217.968495 s. 289 ,Returning DM Positions 44.060167 s. 290 ,Returning DM Positions 52.912735 s. 291 ,Returning DM Positions 38.166918 s. 292 ,Returning DM Positions 36.652352 s. 293 ,Returning DM Positions 40.879841 s. 294 ,Returning DM Positions 55.896002 s. 295 ,Returning DM Positions 55.597522 s. 296 ,Returning DM Positions 46.889938 s. 297 ,Returning DM Positions 40.104389 s. 298 ,Returning DM Positions 41.991030 s. 299 ,Returning DM Positions 44.919619 s. 300 ,Returning DM Positions 46.684174 s. 301 ,Returning DM Positions 37.251948 s. 302 ,Returning DM Positions 50.415917 s. 303 ,Returning DM Positions 45.222695 s. 304 ,Returning DM Positions 48.057070 s. 305 ,Returning DM Positions 43.529459 s. 306 ,Returning DM Positions 53.103523 s. 307 ,Returning DM Positions 39.935632 s. 308 ,Returning DM Positions 36.525156 s. 309 ,Returning DM Positions 56.932001 s. 310 ,Returning DM Positions 42.916055 s. 311 ,Returning DM Positions 45.810113 s. 312 ,Returning DM Positions 35.129134 s. 313 ,Returning DM Positions 59.197894 s. 314 ,Returning DM Positions 53.954292 s. 315 ,Returning DM Positions 49.895464 s. 316 ,Returning DM Positions 41.434467 s. 317 ,Returning DM Positions 50.571852 s. 318 ,Returning DM Positions 56.538154 s. 319 ,Returning DM Positions 53.343257 s. 320 ,Returning DM Positions 47.940166 s. 321 ,Returning DM Positions 51.524338 s. 322 ,Returning DM Positions 41.065753 s. 323 ,Returning DM Positions 38.949559 s. 324 ,Returning DM Positions 246.950756 s. 325 ,Returning DM Positions 59.091646 s. 326 ,Returning DM Positions 43.737118 s. 327 ,Returning DM Positions 162.143904 s. 328 ,Returning DM Positions 58.525069 s. 329 ,Returning DM Positions 69.630337 s. 330 ,Returning DM Positions 597.547492 s. 331 ,Returning DM Positions 39.909442 s. 332 ,Returning DM Positions 37.901819 s. 333 ,Returning DM Positions 42.080470 s. 334 ,Returning DM Positions 72.723154 s. 335 ,Returning DM Positions 48.330019 s. 336 ,Returning DM Positions 47.115901 s. 337 ,Returning DM Positions 57.834247 s. 338 ,Returning DM Positions 43.865574 s. 339 ,Returning DM Positions 39.928522 s. 340 ,Returning DM Positions 41.426796 s. 341 ,Returning DM Positions 44.103147 s. 342 ,Returning DM Positions 50.012167 s. 343 ,Returning DM Positions 50.200857 s. 344 ,Returning DM Positions 54.025968 s. 345 ,Returning DM Positions 38.883153 s. 346 ,Returning DM Positions 40.310398 s. 347 ,Returning DM Positions 43.385725 s. 348 ,Returning DM Positions 66.529093 s. 349 ,Returning DM Positions 41.003021 s. 350 ,Returning DM Positions 50.514502 s. 351 ,Returning DM Positions 41.805768 s. 352 ,Returning DM Positions 51.759195 s. 353 ,Returning DM Positions 47.436666 s. 354 ,Returning DM Positions 41.895516 s. 355 ,Returning DM Positions 49.258050 s. 356 ,Returning DM Positions 85.479489 s. 357 ,Returning DM Positions 41.636994 s. 358 ,Returning DM Positions 58.769939 s. 359 ,Returning DM Positions 55.615698 s. 360 ,Returning DM Positions 50.489549 s. 361 ,Returning DM Positions 53.961739 s. 362 ,Returning DM Positions 51.825651 s. 363 ,Returning DM Positions 56.749763 s. 364 ,Returning DM Positions 48.405143 s. 365 ,Returning DM Positions 45.755516 s. 366 ,Returning DM Positions 47.752856 s. 367 ,Returning DM Positions 49.772810 s. 368 ,Returning DM Positions 45.917100 s. 369 ,Returning DM Positions 42.502296 s. 370 ,Returning DM Positions 73.581183 s. 371 ,Returning DM Positions 202.341788 s. 372 ,Returning DM Positions 67.130182 s. 373 ,Returning DM Positions 42.658697 s. 374 ,Returning DM Positions 58.249956 s. 375 ,Returning DM Positions 57.307707 s. 376 ,Returning DM Positions 47.064054 s. 377 ,Returning DM Positions 65.842666 s. 378 ,Returning DM Positions 43.384042 s. 379 ,Returning DM Positions 56.418867 s. 380 ,Returning DM Positions 42.851543 s. 381 ,Returning DM Positions 133.885670 s. 382 ,Returning DM Positions 40.774373 s. 383 ,Returning DM Positions 43.206785 s. 384 ,Returning DM Positions 40.502748 s. 385 ,Returning DM Positions 40.130596 s. 386 ,Returning DM Positions 47.389595 s. 387 ,Returning DM Positions 38.118967 s. 388 ,Returning DM Positions 51.014190 s. 389 ,Returning DM Positions 47.897410 s. 390 ,Returning DM Positions 50.603735 s. 391 ,Returning DM Positions 114.637599 s. 392 ,Returning DM Positions 147.744329 s. 393 ,Returning DM Positions 60.247224 s. 394 ,Returning DM Positions 55.273425 s. 395 ,Returning DM Positions 38.983756 s. 396 ,Returning DM Positions 38.370264 s. 397 ,Returning DM Positions 43.793579 s. 398 ,Returning DM Positions 48.627239 s. 399 ,Returning DM Positions 43.704439 s. 400 ,Returning DM Positions 51.480573 s. 401 ,Returning DM Positions 45.827287 s. 402 ,Returning DM Positions 38.633766 s. 403 ,Returning DM Positions 43.085098 s. 404 ,Returning DM Positions 39.834360 s. 405 ,Returning DM Positions 41.023363 s. 406 ,Returning DM Positions 41.562496 s. 407 ,Returning DM Positions 49.547254 s. 408 ,Returning DM Positions 47.626842 s. 409 ,Returning DM Positions 48.377432 s. 410 ,Returning DM Positions 53.127023 s. 411 ,Returning DM Positions 49.713989 s. 412 ,Returning DM Positions 51.118086 s. 413 ,Returning DM Positions 43.861520 s. 414 ,Returning DM Positions 44.653101 s. 415 ,Returning DM Positions 55.682187 s. 416 ,Returning DM Positions 45.376061 s. 417 ,Returning DM Positions 42.967053 s. 418 ,Returning DM Positions 55.264381 s. 419 ,Returning DM Positions 47.097006 s. 420 ,Returning DM Positions 46.985958 s. 421 ,Returning DM Positions 45.260533 s. 422 ,Returning DM Positions 53.197648 s. 423 ,Returning DM Positions 55.955639 s. 424 ,Returning DM Positions 61.489531 s. 425 ,Returning DM Positions 72.855227 s. 426 ,Returning DM Positions 52.754157 s. 427 ,Returning DM Positions 62.864767 s. 428 ,Returning DM Positions 54.353206 s. 429 ,Returning DM Positions 51.403277 s. 430 ,Returning DM Positions 43.595654 s. 431 ,Returning DM Positions 54.638009 s. 432 ,Returning DM Positions 49.758974 s. 433 ,Returning DM Positions 44.392011 s. 434 ,Returning DM Positions 58.581956 s. 435 ,Returning DM Positions 46.259884 s. 436 ,Returning DM Positions 50.027617 s. 437 ,Returning DM Positions 581.640156 s. 438 ,Returning DM Positions 49.566267 s. 439 ,Returning DM Positions 46.006141 s. 440 ,Returning DM Positions 46.730331 s. 441 ,Returning DM Positions 36.692934 s. 442 ,Returning DM Positions 46.859017 s. 443 ,Returning DM Positions 45.989913 s. 444 ,Returning DM Positions 39.755156 s. 445 ,Returning DM Positions 62.344973 s. 446 ,Returning DM Positions 192.250703 s. 447 ,Returning DM Positions 44.202347 s. 448 ,Returning DM Positions 53.807785 s. 449 ,Returning DM Positions 47.787614 s. 450 ,Returning DM Positions 69.225250 s. 451 ,Returning DM Positions 64.831177 s. 452 ,Returning DM Positions 41.811014 s. 453 ,Returning DM Positions 43.038856 s. 454 ,Returning DM Positions 44.175777 s. 455 ,Returning DM Positions 92.222251 s. 456 ,Returning DM Positions 89.045178 s. 457 ,Returning DM Positions 44.290426 s. 458 ,Returning DM Positions 42.105298 s. 459 ,Returning DM Positions 48.990762 s. 460 ,Returning DM Positions 47.260827 s. 461 ,Returning DM Positions 46.133823 s. 462 ,Returning DM Positions 47.995849 s. 463 ,Returning DM Positions 47.233761 s. 464 ,Returning DM Positions 43.841338 s. 465 ,Returning DM Positions 48.231867 s. 466 ,Returning DM Positions 43.431142 s. 467 ,Returning DM Positions 122.400445 s. 468 ,Returning DM Positions 44.959718 s. 469 ,Returning DM Positions 47.187552 s. 470 ,Returning DM Positions 58.281118 s. 471 ,Returning DM Positions 54.193462 s. 472 ,Returning DM Positions 63.466470 s. 473 ,Returning DM Positions 47.159656 s. 474 ,Returning DM Positions 43.744367 s. 475 ,Returning DM Positions 44.002954 s. 476 ,Returning DM Positions 70.799842 s. 477 ,Returning DM Positions 56.608969 s. 478 ,Returning DM Positions 46.834882 s. 479 ,Returning DM Positions 46.179053 s. 480 ,Returning DM Positions 45.102349 s. 481 ,Returning DM Positions 59.453735 s. 482 ,Returning DM Positions 45.379437 s. 483 ,Returning DM Positions 54.512213 s. 484 ,Returning DM Positions 43.307754 s. 485 ,Returning DM Positions 132.868501 s. 486 ,Returning DM Positions 195.463878 s. 487 ,Returning DM Positions 44.080508 s. 488 ,Returning DM Positions 48.438794 s. 489 ,Returning DM Positions 44.358165 s. 490 ,Returning DM Positions 57.201629 s. 491 ,Returning DM Positions 50.232809 s. 492 ,Returning DM Positions 41.923515 s. 493 ,Returning DM Positions 50.990171 s. 494 ,Returning DM Positions 50.735100 s. 495 ,Returning DM Positions 42.662746 s. 496 ,Returning DM Positions 55.451264 s. 497 ,Returning DM Positions 60.527736 s. 498 ,Returning DM Positions 159.048244 s. 499 ,Returning DM Positions 41.764618 s. 500 ,Returning DM Positions 47.820993 s. 501 ,Returning DM Positions 47.524187 s. 502 ,Returning DM Positions 44.975470 s. 503 ,Returning DM Positions 60.401641 s. 504 ,Returning DM Positions 62.688514 s. 505 ,Returning DM Positions 51.865276 s. 506 ,Returning DM Positions 44.745231 s. 507 ,Returning DM Positions 54.323832 s. 508 ,Returning DM Positions 40.320305 s. 509 ,Returning DM Positions 55.214127 s. 510 ,Returning DM Positions 38.482095 s. 511 ,Returning DM Positions 42.977619 s.
In [119]:
stack.shape
Out[119]:
(12288,)
In [140]:
np.save(dataDir+'millHealpix2Mpc.npy',stack)
In [120]:
hp.mollview(stack, title="Mollview image RING")
hp.graticule()
0.0 180.0 -180.0 180.0 The interval between parallels is 30 deg -0.00'. The interval between meridians is 30 deg -0.00'.
In [121]:
hp.mollview(hp.sphtfunc.smoothing(stack,0.1), title="Mollview image RING")
hp.graticule()
Sigma is 145.987664 arcmin (0.042466 rad) -> fwhm is 343.774677 arcmin Sigma is 0.000000 arcmin (0.000000 rad) -> fwhm is 0.000000 arcmin 0.0 180.0 -180.0 180.0 The interval between parallels is 30 deg -0.00'. The interval between meridians is 30 deg -0.00'.
Power spec¶
In [122]:
LMAX = 2*NSIDE
cl = hp.anafast(stack, lmax=LMAX)
ell = np.arange(len(cl))
In [123]:
alm = hp.sphtfunc.map2alm(stack, lmax=LMAX)
In [124]:
plt.figure(figsize=(10, 5))
plt.plot(ell, ell * (ell + 1) * cl)
plt.xlabel("$\ell$")
plt.ylabel("$\ell(\ell+1)C_{\ell}$")
plt.grid()
In [126]:
plt.figure(figsize=(10, 5))
k = 17
plt.plot(ell[:k], ell[:k] * (ell[:k] + 1) * cl[:k])
plt.xlabel("$\ell$")
plt.ylabel("$\ell(\ell+1)C_{\ell}$")
plt.grid()
In [127]:
LIM = 10000
plt.plot(alm.real,alm.imag,'.')
plt.xlim(-LIM,LIM)
plt.ylim(-LIM,LIM)
Out[127]:
(-10000, 10000)
In [135]:
plt.semilogy(np.abs(alm.real[:64]))
Out[135]:
[<matplotlib.lines.Line2D at 0x7f81c19faef0>]
In [132]:
plt.semilogy(np.abs(alm.imag[:100]))
Out[132]:
[<matplotlib.lines.Line2D at 0x7f81d80c1048>]
R= 0.5 Mpc/h¶
In [147]:
stackRadius = 0.5
#for i in np.arange(32,dtype=int):
for i in np.arange(480,dtype=int)+32:
print(i,",",end='')
start = time.time()
if i == 0:
stack05 = healPixStat(i)
else:
stack05 += healPixStat(i)
end = time.time()
print("%f s." %(end-start))
32 ,Returning DM Positions 26.640046 s. 33 ,Returning DM Positions 276.800972 s. 34 ,Returning DM Positions 44.910847 s. 35 ,Returning DM Positions 27.823808 s. 36 ,Returning DM Positions 28.407546 s. 37 ,Returning DM Positions 35.103326 s. 38 ,Returning DM Positions 39.755765 s. 39 ,Returning DM Positions 98.826003 s. 40 ,Returning DM Positions 22.404597 s. 41 ,Returning DM Positions 22.648284 s. 42 ,Returning DM Positions 23.031185 s. 43 ,Returning DM Positions 25.616002 s. 44 ,Returning DM Positions 38.712941 s. 45 ,Returning DM Positions 22.395253 s. 46 ,Returning DM Positions 29.846645 s. 47 ,Returning DM Positions 24.937054 s. 48 ,Returning DM Positions 31.408854 s. 49 ,Returning DM Positions 28.921843 s. 50 ,Returning DM Positions 23.813789 s. 51 ,Returning DM Positions 24.052422 s. 52 ,Returning DM Positions 34.556973 s. 53 ,Returning DM Positions 44.688859 s. 54 ,Returning DM Positions 70.955589 s. 55 ,Returning DM Positions 29.642818 s. 56 ,Returning DM Positions 33.932867 s. 57 ,Returning DM Positions 35.233427 s. 58 ,Returning DM Positions 29.573798 s. 59 ,Returning DM Positions 22.935790 s. 60 ,Returning DM Positions 21.807802 s. 61 ,Returning DM Positions 22.219004 s. 62 ,Returning DM Positions 41.102638 s. 63 ,Returning DM Positions 27.918973 s. 64 ,Returning DM Positions 21.762643 s. 65 ,Returning DM Positions 26.980055 s. 66 ,Returning DM Positions 115.885428 s. 67 ,Returning DM Positions 25.802429 s. 68 ,Returning DM Positions 44.724934 s. 69 ,Returning DM Positions 50.333103 s. 70 ,Returning DM Positions 22.173740 s. 71 ,Returning DM Positions 26.158097 s. 72 ,Returning DM Positions 28.767570 s. 73 ,Returning DM Positions 42.468791 s. 74 ,Returning DM Positions 33.517922 s. 75 ,Returning DM Positions 24.603220 s. 76 ,Returning DM Positions 22.526309 s. 77 ,Returning DM Positions 34.950022 s. 78 ,Returning DM Positions 35.725925 s. 79 ,Returning DM Positions 26.978832 s. 80 ,Returning DM Positions 42.742816 s. 81 ,Returning DM Positions 43.620402 s. 82 ,Returning DM Positions 41.071480 s. 83 ,Returning DM Positions 22.978714 s. 84 ,Returning DM Positions 26.701724 s. 85 ,Returning DM Positions 28.082726 s. 86 ,Returning DM Positions 23.264729 s. 87 ,Returning DM Positions 23.899321 s. 88 ,Returning DM Positions 30.832002 s. 89 ,Returning DM Positions 30.952905 s. 90 ,Returning DM Positions 23.424645 s. 91 ,Returning DM Positions 33.848742 s. 92 ,Returning DM Positions 36.013268 s. 93 ,Returning DM Positions 33.157206 s. 94 ,Returning DM Positions 23.041985 s. 95 ,Returning DM Positions 23.640386 s. 96 ,Returning DM Positions 28.673898 s. 97 ,Returning DM Positions 41.065223 s. 98 ,Returning DM Positions 26.669877 s. 99 ,Returning DM Positions 44.401785 s. 100 ,Returning DM Positions 35.689551 s. 101 ,Returning DM Positions 37.013549 s. 102 ,Returning DM Positions 250.930023 s. 103 ,Returning DM Positions 28.529711 s. 104 ,Returning DM Positions 24.680299 s. 105 ,Returning DM Positions 23.386260 s. 106 ,Returning DM Positions 23.897731 s. 107 ,Returning DM Positions 25.606868 s. 108 ,Returning DM Positions 26.201368 s. 109 ,Returning DM Positions 24.668138 s. 110 ,Returning DM Positions 26.388054 s. 111 ,Returning DM Positions 31.666105 s. 112 ,Returning DM Positions 38.034609 s. 113 ,Returning DM Positions 34.412947 s. 114 ,Returning DM Positions 24.499499 s. 115 ,Returning DM Positions 30.954450 s. 116 ,Returning DM Positions 39.881712 s. 117 ,Returning DM Positions 37.065298 s. 118 ,Returning DM Positions 33.070819 s. 119 ,Returning DM Positions 27.300228 s. 120 ,Returning DM Positions 32.036485 s. 121 ,Returning DM Positions 31.059719 s. 122 ,Returning DM Positions 27.113439 s. 123 ,Returning DM Positions 23.595489 s. 124 ,Returning DM Positions 32.106075 s. 125 ,Returning DM Positions 25.406493 s. 126 ,Returning DM Positions 27.822399 s. 127 ,Returning DM Positions 27.422545 s. 128 ,Returning DM Positions 23.764619 s. 129 ,Returning DM Positions 25.910247 s. 130 ,Returning DM Positions 45.269802 s. 131 ,Returning DM Positions 29.602110 s. 132 ,Returning DM Positions 22.883723 s. 133 ,Returning DM Positions 42.938592 s. 134 ,Returning DM Positions 43.647199 s. 135 ,Returning DM Positions 39.526592 s. 136 ,Returning DM Positions 37.330332 s. 137 ,Returning DM Positions 30.837323 s. 138 ,Returning DM Positions 785.901667 s. 139 ,Returning DM Positions 22.892281 s. 140 ,Returning DM Positions 23.149672 s. 141 ,Returning DM Positions 32.656651 s. 142 ,Returning DM Positions 778.357016 s. 143 ,Returning DM Positions 44.163966 s. 144 ,Returning DM Positions 33.707136 s. 145 ,Returning DM Positions 41.110011 s. 146 ,Returning DM Positions 35.477899 s. 147 ,Returning DM Positions 24.207496 s. 148 ,Returning DM Positions 22.801279 s. 149 ,Returning DM Positions 25.851425 s. 150 ,Returning DM Positions 26.679110 s. 151 ,Returning DM Positions 22.471761 s. 152 ,Returning DM Positions 22.144336 s. 153 ,Returning DM Positions 30.278884 s. 154 ,Returning DM Positions 33.081502 s. 155 ,Returning DM Positions 27.869640 s. 156 ,Returning DM Positions 43.107902 s. 157 ,Returning DM Positions 110.627278 s. 158 ,Returning DM Positions 39.960916 s. 159 ,Returning DM Positions 23.091312 s. 160 ,Returning DM Positions 22.412832 s. 161 ,Returning DM Positions 40.670534 s. 162 ,Returning DM Positions 28.394225 s. 163 ,Returning DM Positions 32.696434 s. 164 ,Returning DM Positions 30.216226 s. 165 ,Returning DM Positions 30.355738 s. 166 ,Returning DM Positions 25.894771 s. 167 ,Returning DM Positions 23.900803 s. 168 ,Returning DM Positions 22.736587 s. 169 ,Returning DM Positions 23.108326 s. 170 ,Returning DM Positions 28.528470 s. 171 ,Returning DM Positions 28.470143 s. 172 ,Returning DM Positions 24.707910 s. 173 ,Returning DM Positions 24.956751 s. 174 ,Returning DM Positions 23.587041 s. 175 ,Returning DM Positions 156.134762 s. 176 ,Returning DM Positions 26.233643 s. 177 ,Returning DM Positions 23.423659 s. 178 ,Returning DM Positions 22.343372 s. 179 ,Returning DM Positions 40.133934 s. 180 ,Returning DM Positions 27.336888 s. 181 ,Returning DM Positions 35.523607 s. 182 ,Returning DM Positions 30.387537 s. 183 ,Returning DM Positions 22.254230 s. 184 ,Returning DM Positions 29.672786 s. 185 ,Returning DM Positions 30.521029 s. 186 ,Returning DM Positions 26.766447 s. 187 ,Returning DM Positions 22.780081 s. 188 ,Returning DM Positions 30.725575 s. 189 ,Returning DM Positions 24.374216 s. 190 ,Returning DM Positions 42.107317 s. 191 ,Returning DM Positions 21.628220 s. 192 ,Returning DM Positions 29.001789 s. 193 ,Returning DM Positions 23.600015 s. 194 ,Returning DM Positions 27.794899 s. 195 ,Returning DM Positions 27.355396 s. 196 ,Returning DM Positions 38.185125 s. 197 ,Returning DM Positions 40.795144 s. 198 ,Returning DM Positions 29.736820 s. 199 ,Returning DM Positions 210.069188 s. 200 ,Returning DM Positions 29.374458 s. 201 ,Returning DM Positions 40.104386 s. 202 ,Returning DM Positions 29.200665 s. 203 ,Returning DM Positions 23.876608 s. 204 ,Returning DM Positions 27.246826 s. 205 ,Returning DM Positions 43.625709 s. 206 ,Returning DM Positions 27.524782 s. 207 ,Returning DM Positions 23.456733 s. 208 ,Returning DM Positions 27.508219 s. 209 ,Returning DM Positions 38.604399 s. 210 ,Returning DM Positions 153.161551 s. 211 ,Returning DM Positions 31.704715 s. 212 ,Returning DM Positions 28.215551 s. 213 ,Returning DM Positions 29.198565 s. 214 ,Returning DM Positions 21.427918 s. 215 ,Returning DM Positions 42.769379 s. 216 ,Returning DM Positions 44.012275 s. 217 ,Returning DM Positions 46.140920 s. 218 ,Returning DM Positions 34.562119 s. 219 ,Returning DM Positions 22.724707 s. 220 ,Returning DM Positions 28.416310 s. 221 ,Returning DM Positions 25.188559 s. 222 ,Returning DM Positions 308.143253 s. 223 ,Returning DM Positions 23.434732 s. 224 ,Returning DM Positions 35.055851 s. 225 ,Returning DM Positions 24.263458 s. 226 ,Returning DM Positions 33.140891 s. 227 ,Returning DM Positions 26.721622 s. 228 ,Returning DM Positions 25.790563 s. 229 ,Returning DM Positions 25.983573 s. 230 ,Returning DM Positions 44.607263 s. 231 ,Returning DM Positions 579.833684 s. 232 ,Returning DM Positions 24.731103 s. 233 ,Returning DM Positions 21.762820 s. 234 ,Returning DM Positions 25.489047 s. 235 ,Returning DM Positions 21.512637 s. 236 ,Returning DM Positions 28.108389 s. 237 ,Returning DM Positions 23.523613 s. 238 ,Returning DM Positions 35.811360 s. 239 ,Returning DM Positions 43.827294 s. 240 ,Returning DM Positions 30.065436 s. 241 ,Returning DM Positions 40.729738 s. 242 ,Returning DM Positions 39.224017 s. 243 ,Returning DM Positions 83.702603 s. 244 ,Returning DM Positions 29.759614 s. 245 ,Returning DM Positions 24.126797 s. 246 ,Returning DM Positions 65.705065 s. 247 ,Returning DM Positions 37.609440 s. 248 ,Returning DM Positions 24.392134 s. 249 ,Returning DM Positions 113.960801 s. 250 ,Returning DM Positions 30.999708 s. 251 ,Returning DM Positions 43.482085 s. 252 ,Returning DM Positions 22.863635 s. 253 ,Returning DM Positions 29.635377 s. 254 ,Returning DM Positions 22.345612 s. 255 ,Returning DM Positions 248.495396 s. 256 ,Returning DM Positions 32.288474 s. 257 ,Returning DM Positions 30.623656 s. 258 ,Returning DM Positions 28.903864 s. 259 ,Returning DM Positions 55.365592 s. 260 ,Returning DM Positions 32.279464 s. 261 ,Returning DM Positions 23.019331 s. 262 ,Returning DM Positions 24.285660 s. 263 ,Returning DM Positions 27.824371 s. 264 ,Returning DM Positions 314.055044 s. 265 ,Returning DM Positions 39.136231 s. 266 ,Returning DM Positions 22.685656 s. 267 ,Returning DM Positions 25.610482 s. 268 ,Returning DM Positions 23.937661 s. 269 ,Returning DM Positions 28.111798 s. 270 ,Returning DM Positions 44.593941 s. 271 ,Returning DM Positions 33.996581 s. 272 ,Returning DM Positions 27.341901 s. 273 ,Returning DM Positions 37.101196 s. 274 ,Returning DM Positions 23.795679 s. 275 ,Returning DM Positions 31.408319 s. 276 ,Returning DM Positions 23.738083 s. 277 ,Returning DM Positions 31.576870 s. 278 ,Returning DM Positions 44.348366 s. 279 ,Returning DM Positions 325.179186 s. 280 ,Returning DM Positions 28.463218 s. 281 ,Returning DM Positions 34.210695 s. 282 ,Returning DM Positions 24.690562 s. 283 ,Returning DM Positions 73.665583 s. 284 ,Returning DM Positions 31.901395 s. 285 ,Returning DM Positions 29.629747 s. 286 ,Returning DM Positions 25.450993 s. 287 ,Returning DM Positions 34.024754 s. 288 ,Returning DM Positions 200.883642 s. 289 ,Returning DM Positions 23.526726 s. 290 ,Returning DM Positions 27.636492 s. 291 ,Returning DM Positions 23.355467 s. 292 ,Returning DM Positions 20.809435 s. 293 ,Returning DM Positions 23.670634 s. 294 ,Returning DM Positions 36.255004 s. 295 ,Returning DM Positions 35.807983 s. 296 ,Returning DM Positions 28.617650 s. 297 ,Returning DM Positions 24.291156 s. 298 ,Returning DM Positions 23.979865 s. 299 ,Returning DM Positions 23.873365 s. 300 ,Returning DM Positions 27.386177 s. 301 ,Returning DM Positions 22.553249 s. 302 ,Returning DM Positions 35.090518 s. 303 ,Returning DM Positions 23.924215 s. 304 ,Returning DM Positions 33.815835 s. 305 ,Returning DM Positions 27.392072 s. 306 ,Returning DM Positions 40.067554 s. 307 ,Returning DM Positions 23.222255 s. 308 ,Returning DM Positions 21.686761 s. 309 ,Returning DM Positions 41.235573 s. 310 ,Returning DM Positions 23.637626 s. 311 ,Returning DM Positions 23.025406 s. 312 ,Returning DM Positions 22.429263 s. 313 ,Returning DM Positions 36.203162 s. 314 ,Returning DM Positions 33.516577 s. 315 ,Returning DM Positions 26.351320 s. 316 ,Returning DM Positions 23.018684 s. 317 ,Returning DM Positions 32.312916 s. 318 ,Returning DM Positions 41.164274 s. 319 ,Returning DM Positions 37.002037 s. 320 ,Returning DM Positions 32.050576 s. 321 ,Returning DM Positions 33.343778 s. 322 ,Returning DM Positions 25.231333 s. 323 ,Returning DM Positions 23.331003 s. 324 ,Returning DM Positions 228.947153 s. 325 ,Returning DM Positions 40.675937 s. 326 ,Returning DM Positions 25.145474 s. 327 ,Returning DM Positions 148.642315 s. 328 ,Returning DM Positions 40.550559 s. 329 ,Returning DM Positions 47.381403 s. 330 ,Returning DM Positions 586.333717 s. 331 ,Returning DM Positions 21.304919 s. 332 ,Returning DM Positions 21.973084 s. 333 ,Returning DM Positions 26.356428 s. 334 ,Returning DM Positions 56.050798 s. 335 ,Returning DM Positions 27.746279 s. 336 ,Returning DM Positions 23.955149 s. 337 ,Returning DM Positions 38.605053 s. 338 ,Returning DM Positions 24.457052 s. 339 ,Returning DM Positions 21.891992 s. 340 ,Returning DM Positions 22.802041 s. 341 ,Returning DM Positions 25.120384 s. 342 ,Returning DM Positions 30.005736 s. 343 ,Returning DM Positions 28.834265 s. 344 ,Returning DM Positions 37.200743 s. 345 ,Returning DM Positions 23.143564 s. 346 ,Returning DM Positions 22.478575 s. 347 ,Returning DM Positions 23.941057 s. 348 ,Returning DM Positions 44.006586 s. 349 ,Returning DM Positions 22.651785 s. 350 ,Returning DM Positions 30.908654 s. 351 ,Returning DM Positions 21.937927 s. 352 ,Returning DM Positions 35.714013 s. 353 ,Returning DM Positions 23.057574 s. 354 ,Returning DM Positions 21.233676 s. 355 ,Returning DM Positions 27.107290 s. 356 ,Returning DM Positions 68.140269 s. 357 ,Returning DM Positions 23.038931 s. 358 ,Returning DM Positions 43.131490 s. 359 ,Returning DM Positions 30.301029 s. 360 ,Returning DM Positions 29.140794 s. 361 ,Returning DM Positions 35.207597 s. 362 ,Returning DM Positions 34.130757 s. 363 ,Returning DM Positions 30.164005 s. 364 ,Returning DM Positions 26.074113 s. 365 ,Returning DM Positions 25.819337 s. 366 ,Returning DM Positions 33.719177 s. 367 ,Returning DM Positions 30.139164 s. 368 ,Returning DM Positions 24.739382 s. 369 ,Returning DM Positions 22.449877 s. 370 ,Returning DM Positions 53.138590 s. 371 ,Returning DM Positions 170.405916 s. 372 ,Returning DM Positions 47.144746 s. 373 ,Returning DM Positions 25.907318 s. 374 ,Returning DM Positions 41.774745 s. 375 ,Returning DM Positions 40.078874 s. 376 ,Returning DM Positions 28.236837 s. 377 ,Returning DM Positions 45.896935 s. 378 ,Returning DM Positions 23.366199 s. 379 ,Returning DM Positions 36.916776 s. 380 ,Returning DM Positions 22.958097 s. 381 ,Returning DM Positions 114.686608 s. 382 ,Returning DM Positions 23.430743 s. 383 ,Returning DM Positions 25.999882 s. 384 ,Returning DM Positions 23.617490 s. 385 ,Returning DM Positions 21.611762 s. 386 ,Returning DM Positions 25.822118 s. 387 ,Returning DM Positions 23.913068 s. 388 ,Returning DM Positions 22.516476 s. 389 ,Returning DM Positions 26.682435 s. 390 ,Returning DM Positions 26.013273 s. 391 ,Returning DM Positions 100.232348 s. 392 ,Returning DM Positions 129.307714 s. 393 ,Returning DM Positions 44.997933 s. 394 ,Returning DM Positions 29.492798 s. 395 ,Returning DM Positions 22.381776 s. 396 ,Returning DM Positions 22.796034 s. 397 ,Returning DM Positions 28.993333 s. 398 ,Returning DM Positions 33.790703 s. 399 ,Returning DM Positions 23.048437 s. 400 ,Returning DM Positions 30.292953 s. 401 ,Returning DM Positions 29.217277 s. 402 ,Returning DM Positions 22.669274 s. 403 ,Returning DM Positions 26.812469 s. 404 ,Returning DM Positions 23.696702 s. 405 ,Returning DM Positions 26.042211 s. 406 ,Returning DM Positions 23.115365 s. 407 ,Returning DM Positions 27.712800 s. 408 ,Returning DM Positions 27.954760 s. 409 ,Returning DM Positions 28.057550 s. 410 ,Returning DM Positions 38.601173 s. 411 ,Returning DM Positions 31.920707 s. 412 ,Returning DM Positions 27.908686 s. 413 ,Returning DM Positions 22.721758 s. 414 ,Returning DM Positions 24.342788 s. 415 ,Returning DM Positions 39.864808 s. 416 ,Returning DM Positions 29.555138 s. 417 ,Returning DM Positions 24.749380 s. 418 ,Returning DM Positions 32.954510 s. 419 ,Returning DM Positions 30.490855 s. 420 ,Returning DM Positions 28.621007 s. 421 ,Returning DM Positions 25.973789 s. 422 ,Returning DM Positions 37.380746 s. 423 ,Returning DM Positions 38.357228 s. 424 ,Returning DM Positions 44.726959 s. 425 ,Returning DM Positions 44.811968 s. 426 ,Returning DM Positions 25.607250 s. 427 ,Returning DM Positions 43.503640 s. 428 ,Returning DM Positions 34.180938 s. 429 ,Returning DM Positions 25.708098 s. 430 ,Returning DM Positions 23.103840 s. 431 ,Returning DM Positions 25.146101 s. 432 ,Returning DM Positions 28.096846 s. 433 ,Returning DM Positions 25.783873 s. 434 ,Returning DM Positions 41.113972 s. 435 ,Returning DM Positions 25.632344 s. 436 ,Returning DM Positions 24.290925 s. 437 ,Returning DM Positions 550.445289 s. 438 ,Returning DM Positions 26.215366 s. 439 ,Returning DM Positions 27.182788 s. 440 ,Returning DM Positions 28.589456 s. 441 ,Returning DM Positions 21.805262 s. 442 ,Returning DM Positions 29.933153 s. 443 ,Returning DM Positions 24.892926 s. 444 ,Returning DM Positions 23.289715 s. 445 ,Returning DM Positions 48.290253 s. 446 ,Returning DM Positions 215.934823 s. 447 ,Returning DM Positions 29.477761 s. 448 ,Returning DM Positions 30.334658 s. 449 ,Returning DM Positions 25.683424 s. 450 ,Returning DM Positions 53.301678 s. 451 ,Returning DM Positions 53.200514 s. 452 ,Returning DM Positions 25.793684 s. 453 ,Returning DM Positions 24.869054 s. 454 ,Returning DM Positions 25.895241 s. 455 ,Returning DM Positions 67.999064 s. 456 ,Returning DM Positions 75.600666 s. 457 ,Returning DM Positions 23.637960 s. 458 ,Returning DM Positions 24.762969 s. 459 ,Returning DM Positions 24.402685 s. 460 ,Returning DM Positions 28.714165 s. 461 ,Returning DM Positions 26.901895 s. 462 ,Returning DM Positions 30.474183 s. 463 ,Returning DM Positions 31.125244 s. 464 ,Returning DM Positions 26.639100 s. 465 ,Returning DM Positions 26.729225 s. 466 ,Returning DM Positions 25.574424 s. 467 ,Returning DM Positions 124.538229 s. 468 ,Returning DM Positions 26.968705 s. 469 ,Returning DM Positions 23.909688 s. 470 ,Returning DM Positions 37.039260 s. 471 ,Returning DM Positions 31.695505 s. 472 ,Returning DM Positions 45.146810 s. 473 ,Returning DM Positions 25.526428 s. 474 ,Returning DM Positions 26.461315 s. 475 ,Returning DM Positions 28.432795 s. 476 ,Returning DM Positions 56.836948 s. 477 ,Returning DM Positions 39.044946 s. 478 ,Returning DM Positions 29.121591 s. 479 ,Returning DM Positions 23.673788 s. 480 ,Returning DM Positions 23.552225 s. 481 ,Returning DM Positions 36.936019 s. 482 ,Returning DM Positions 26.241528 s. 483 ,Returning DM Positions 34.337635 s. 484 ,Returning DM Positions 24.636920 s. 485 ,Returning DM Positions 124.826536 s. 486 ,Returning DM Positions 209.252087 s. 487 ,Returning DM Positions 25.832546 s. 488 ,Returning DM Positions 31.087046 s. 489 ,Returning DM Positions 26.303778 s. 490 ,Returning DM Positions 39.995721 s. 491 ,Returning DM Positions 33.648493 s. 492 ,Returning DM Positions 27.847062 s. 493 ,Returning DM Positions 37.768825 s. 494 ,Returning DM Positions 39.095352 s. 495 ,Returning DM Positions 25.066275 s. 496 ,Returning DM Positions 35.006194 s. 497 ,Returning DM Positions 36.431961 s. 498 ,Returning DM Positions 144.579519 s. 499 ,Returning DM Positions 23.252688 s. 500 ,Returning DM Positions 25.243801 s. 501 ,Returning DM Positions 25.447943 s. 502 ,Returning DM Positions 21.851802 s. 503 ,Returning DM Positions 37.507684 s. 504 ,Returning DM Positions 41.621125 s. 505 ,Returning DM Positions 37.370099 s. 506 ,Returning DM Positions 28.519063 s. 507 ,Returning DM Positions 39.412805 s. 508 ,Returning DM Positions 23.642006 s. 509 ,Returning DM Positions 38.149780 s. 510 ,Returning DM Positions 22.462818 s. 511 ,Returning DM Positions 25.406737 s.
In [148]:
hp.mollview(stack05, title="Mollview image RING")
hp.graticule()
0.0 180.0 -180.0 180.0 The interval between parallels is 30 deg -0.00'. The interval between meridians is 30 deg -0.00'.
In [149]:
LMAX = 2*NSIDE
cl05 = hp.anafast(stack05, lmax=LMAX)
ell05 = np.arange(len(cl05))
In [154]:
plt.figure(figsize=(10, 5))
plt.plot(ell05, ell05 * (ell05 + 1) * cl05)
plt.plot(ell, ell * (ell + 1) * cl/2000)
plt.xlabel("$\ell$")
plt.ylabel("$\ell(\ell+1)C_{\ell}$")
plt.title('R = 0.5 Mpc/h')
plt.grid()
#plt.xlim(0,17)
R= 4 x halfMassRadius¶
In [158]:
stackRadius = 0.5
#for i in np.arange(5,dtype=int):
for i in np.arange(512,dtype=int):
print(i,",",end='')
start = time.time()
if i == 0:
stackR = healPixStatHalfMassR(i)
else:
stackR += healPixStatHalfMassR(i)
end = time.time()
print("%f s." %(end-start))
0 ,Returning DM Positions 24.575250 s. 1 ,Returning DM Positions 23.364319 s. 2 ,Returning DM Positions 27.091832 s. 3 ,Returning DM Positions 41.781519 s. 4 ,Returning DM Positions 39.324396 s. 5 ,Returning DM Positions 23.255033 s. 6 ,Returning DM Positions 26.803133 s. 7 ,Returning DM Positions 187.491047 s. 8 ,Returning DM Positions 29.692605 s. 9 ,Returning DM Positions 43.679567 s. 10 ,Returning DM Positions 26.063485 s. 11 ,Returning DM Positions 26.294942 s. 12 ,Returning DM Positions 30.931000 s. 13 ,Returning DM Positions 21.325915 s. 14 ,Returning DM Positions 27.236347 s. 15 ,Returning DM Positions 30.910559 s. 16 ,Returning DM Positions 27.051880 s. 17 ,Returning DM Positions 36.236677 s. 18 ,Returning DM Positions 37.931621 s. 19 ,Returning DM Positions 30.885365 s. 20 ,Returning DM Positions 37.113486 s. 21 ,Returning DM Positions 35.873267 s. 22 ,Returning DM Positions 29.914214 s. 23 ,Returning DM Positions 26.695277 s. 24 ,Returning DM Positions 30.631353 s. 25 ,Returning DM Positions 24.485146 s. 26 ,Returning DM Positions 32.843754 s. 27 ,Returning DM Positions 23.459400 s. 28 ,Returning DM Positions 22.493772 s. 29 ,Returning DM Positions 35.305626 s. 30 ,Returning DM Positions 34.889096 s. 31 ,Returning DM Positions 28.160085 s. 32 ,Returning DM Positions 26.766532 s. 33 ,Returning DM Positions 282.174401 s. 34 ,Returning DM Positions 49.484938 s. 35 ,Returning DM Positions 28.217067 s. 36 ,Returning DM Positions 29.858187 s. 37 ,Returning DM Positions 34.667774 s. 38 ,Returning DM Positions 40.539066 s. 39 ,Returning DM Positions 99.458121 s. 40 ,Returning DM Positions 22.647775 s. 41 ,Returning DM Positions 22.496357 s. 42 ,Returning DM Positions 22.725809 s. 43 ,Returning DM Positions 24.539719 s. 44 ,Returning DM Positions 37.956964 s. 45 ,Returning DM Positions 22.314800 s. 46 ,Returning DM Positions 30.845907 s. 47 ,Returning DM Positions 25.071298 s. 48 ,Returning DM Positions 32.063923 s. 49 ,Returning DM Positions 28.185606 s. 50 ,Returning DM Positions 24.067510 s. 51 ,Returning DM Positions 24.600673 s. 52 ,Returning DM Positions 36.582580 s. 53 ,Returning DM Positions 44.965494 s. 54 ,Returning DM Positions 70.411319 s. 55 ,Returning DM Positions 29.911511 s. 56 ,Returning DM Positions 34.678796 s. 57 ,Returning DM Positions 35.303141 s. 58 ,Returning DM Positions 29.719052 s. 59 ,Returning DM Positions 23.236583 s. 60 ,Returning DM Positions 22.411567 s. 61 ,Returning DM Positions 21.896923 s. 62 ,Returning DM Positions 39.265411 s. 63 ,Returning DM Positions 27.475785 s. 64 ,Returning DM Positions 21.597541 s. 65 ,Returning DM Positions 27.215485 s. 66 ,Returning DM Positions 117.204368 s. 67 ,Returning DM Positions 25.581069 s. 68 ,Returning DM Positions 46.815684 s. 69 ,Returning DM Positions 50.898006 s. 70 ,Returning DM Positions 22.515559 s. 71 ,Returning DM Positions 25.859993 s. 72 ,Returning DM Positions 29.100125 s. 73 ,Returning DM Positions 42.701078 s. 74 ,Returning DM Positions 33.302404 s. 75 ,Returning DM Positions 25.112358 s. 76 ,Returning DM Positions 21.788934 s. 77 ,Returning DM Positions 33.646089 s. 78 ,Returning DM Positions 34.801901 s. 79 ,Returning DM Positions 24.820335 s. 80 ,Returning DM Positions 43.470672 s. 81 ,Returning DM Positions 43.802568 s. 82 ,Returning DM Positions 41.230489 s. 83 ,Returning DM Positions 23.326346 s. 84 ,Returning DM Positions 27.011699 s. 85 ,Returning DM Positions 28.107424 s. 86 ,Returning DM Positions 24.435101 s. 87 ,Returning DM Positions 24.413422 s. 88 ,Returning DM Positions 30.718047 s. 89 ,Returning DM Positions 32.865079 s. 90 ,Returning DM Positions 23.769451 s. 91 ,Returning DM Positions 34.600225 s. 92 ,Returning DM Positions 36.514998 s. 93 ,Returning DM Positions 33.186440 s. 94 ,Returning DM Positions 23.002081 s. 95 ,Returning DM Positions 22.825885 s. 96 ,Returning DM Positions 30.275384 s. 97 ,Returning DM Positions 41.741939 s. 98 ,Returning DM Positions 24.661641 s. 99 ,Returning DM Positions 45.111823 s. 100 ,Returning DM Positions 36.547996 s. 101 ,Returning DM Positions 38.219899 s. 102 ,Returning DM Positions 256.047485 s. 103 ,Returning DM Positions 31.610405 s. 104 ,Returning DM Positions 28.919739 s. 105 ,Returning DM Positions 24.558046 s. 106 ,Returning DM Positions 24.976120 s. 107 ,Returning DM Positions 26.284880 s. 108 ,Returning DM Positions 27.441189 s. 109 ,Returning DM Positions 25.174846 s. 110 ,Returning DM Positions 27.563340 s. 111 ,Returning DM Positions 30.495456 s. 112 ,Returning DM Positions 38.126711 s. 113 ,Returning DM Positions 34.982844 s. 114 ,Returning DM Positions 25.926380 s. 115 ,Returning DM Positions 30.228293 s. 116 ,Returning DM Positions 39.372534 s. 117 ,Returning DM Positions 37.738398 s. 118 ,Returning DM Positions 34.573642 s. 119 ,Returning DM Positions 27.663980 s. 120 ,Returning DM Positions 31.990406 s. 121 ,Returning DM Positions 31.136579 s. 122 ,Returning DM Positions 27.575219 s. 123 ,Returning DM Positions 23.939588 s. 124 ,Returning DM Positions 33.146774 s. 125 ,Returning DM Positions 25.269595 s. 126 ,Returning DM Positions 28.391950 s. 127 ,Returning DM Positions 27.368140 s. 128 ,Returning DM Positions 24.176176 s. 129 ,Returning DM Positions 26.085346 s. 130 ,Returning DM Positions 44.244957 s. 131 ,Returning DM Positions 28.026293 s. 132 ,Returning DM Positions 23.168900 s. 133 ,Returning DM Positions 42.748762 s. 134 ,Returning DM Positions 43.206405 s. 135 ,Returning DM Positions 39.524085 s. 136 ,Returning DM Positions 39.064190 s. 137 ,Returning DM Positions 32.767586 s. 138 ,Returning DM Positions 807.460489 s. 139 ,Returning DM Positions 23.932396 s. 140 ,Returning DM Positions 24.381042 s. 141 ,Returning DM Positions 31.644276 s. 142 ,Returning DM Positions 803.536655 s. 143 ,Returning DM Positions 44.069793 s. 144 ,Returning DM Positions 34.099026 s. 145 ,Returning DM Positions 43.198670 s. 146 ,Returning DM Positions 35.384427 s. 147 ,Returning DM Positions 23.784215 s. 148 ,Returning DM Positions 22.526107 s. 149 ,Returning DM Positions 25.444463 s. 150 ,Returning DM Positions 26.561212 s. 151 ,Returning DM Positions 22.376010 s. 152 ,Returning DM Positions 21.830090 s. 153 ,Returning DM Positions 31.216666 s. 154 ,Returning DM Positions 33.623892 s. 155 ,Returning DM Positions 27.458173 s. 156 ,Returning DM Positions 43.174638 s. 157 ,Returning DM Positions 106.516704 s. 158 ,Returning DM Positions 40.007625 s. 159 ,Returning DM Positions 22.799504 s. 160 ,Returning DM Positions 21.280553 s. 161 ,Returning DM Positions 39.895397 s. 162 ,Returning DM Positions 30.042650 s. 163 ,Returning DM Positions 33.351804 s. 164 ,Returning DM Positions 30.271392 s. 165 ,Returning DM Positions 29.950456 s. 166 ,Returning DM Positions 25.411694 s. 167 ,Returning DM Positions 23.899385 s. 168 ,Returning DM Positions 23.554075 s. 169 ,Returning DM Positions 22.787381 s. 170 ,Returning DM Positions 29.010229 s. 171 ,Returning DM Positions 28.260373 s. 172 ,Returning DM Positions 24.502982 s. 173 ,Returning DM Positions 24.996908 s. 174 ,Returning DM Positions 22.675084 s. 175 ,Returning DM Positions 149.361643 s. 176 ,Returning DM Positions 26.106400 s. 177 ,Returning DM Positions 22.671659 s. 178 ,Returning DM Positions 22.266970 s. 179 ,Returning DM Positions 39.791958 s. 180 ,Returning DM Positions 28.871027 s. 181 ,Returning DM Positions 36.190286 s. 182 ,Returning DM Positions 30.489052 s. 183 ,Returning DM Positions 21.689847 s. 184 ,Returning DM Positions 29.550951 s. 185 ,Returning DM Positions 30.151718 s. 186 ,Returning DM Positions 26.559699 s. 187 ,Returning DM Positions 22.830055 s. 188 ,Returning DM Positions 30.174374 s. 189 ,Returning DM Positions 24.397134 s. 190 ,Returning DM Positions 41.445339 s. 191 ,Returning DM Positions 21.531821 s. 192 ,Returning DM Positions 28.940079 s. 193 ,Returning DM Positions 23.055987 s. 194 ,Returning DM Positions 27.493493 s. 195 ,Returning DM Positions 25.909903 s. 196 ,Returning DM Positions 36.979729 s. 197 ,Returning DM Positions 39.655869 s. 198 ,Returning DM Positions 29.099476 s. 199 ,Returning DM Positions 221.491353 s. 200 ,Returning DM Positions 28.572854 s. 201 ,Returning DM Positions 39.454130 s. 202 ,Returning DM Positions 28.789186 s. 203 ,Returning DM Positions 23.217117 s. 204 ,Returning DM Positions 27.547004 s. 205 ,Returning DM Positions 43.953397 s. 206 ,Returning DM Positions 27.060159 s. 207 ,Returning DM Positions 23.188770 s. 208 ,Returning DM Positions 25.965309 s. 209 ,Returning DM Positions 37.676711 s. 210 ,Returning DM Positions 153.161816 s. 211 ,Returning DM Positions 32.427165 s. 212 ,Returning DM Positions 27.164065 s. 213 ,Returning DM Positions 29.210632 s. 214 ,Returning DM Positions 22.271254 s. 215 ,Returning DM Positions 42.937204 s. 216 ,Returning DM Positions 42.484279 s. 217 ,Returning DM Positions 45.918970 s. 218 ,Returning DM Positions 34.331460 s. 219 ,Returning DM Positions 22.894716 s. 220 ,Returning DM Positions 28.745949 s. 221 ,Returning DM Positions 25.066995 s. 222 ,Returning DM Positions 307.471565 s. 223 ,Returning DM Positions 23.171467 s. 224 ,Returning DM Positions 35.667512 s. 225 ,Returning DM Positions 24.319431 s. 226 ,Returning DM Positions 32.447159 s. 227 ,Returning DM Positions 25.459307 s. 228 ,Returning DM Positions 26.786920 s. 229 ,Returning DM Positions 26.342339 s. 230 ,Returning DM Positions 44.543947 s. 231 ,Returning DM Positions 598.184579 s. 232 ,Returning DM Positions 24.715746 s. 233 ,Returning DM Positions 21.865848 s. 234 ,Returning DM Positions 25.139390 s. 235 ,Returning DM Positions 21.499960 s. 236 ,Returning DM Positions 26.730746 s. 237 ,Returning DM Positions 22.183360 s. 238 ,Returning DM Positions 35.791458 s. 239 ,Returning DM Positions 43.722462 s. 240 ,Returning DM Positions 29.529824 s. 241 ,Returning DM Positions 43.785486 s. 242 ,Returning DM Positions 38.992019 s. 243 ,Returning DM Positions 83.702572 s. 244 ,Returning DM Positions 29.024495 s. 245 ,Returning DM Positions 24.083737 s. 246 ,Returning DM Positions 64.637150 s. 247 ,Returning DM Positions 37.843870 s. 248 ,Returning DM Positions 24.454255 s. 249 ,Returning DM Positions 112.127660 s. 250 ,Returning DM Positions 30.256332 s. 251 ,Returning DM Positions 42.687321 s. 252 ,Returning DM Positions 22.305283 s. 253 ,Returning DM Positions 29.627974 s. 254 ,Returning DM Positions 21.974687 s. 255 ,Returning DM Positions 252.472817 s. 256 ,Returning DM Positions 31.822058 s. 257 ,Returning DM Positions 30.478046 s. 258 ,Returning DM Positions 27.748078 s. 259 ,Returning DM Positions 55.034884 s. 260 ,Returning DM Positions 33.181679 s. 261 ,Returning DM Positions 21.874576 s. 262 ,Returning DM Positions 23.789087 s. 263 ,Returning DM Positions 26.858021 s. 264 ,Returning DM Positions 316.199030 s. 265 ,Returning DM Positions 38.517058 s. 266 ,Returning DM Positions 22.771643 s. 267 ,Returning DM Positions 25.386533 s. 268 ,Returning DM Positions 23.811099 s. 269 ,Returning DM Positions 27.939110 s. 270 ,Returning DM Positions 44.003875 s. 271 ,Returning DM Positions 33.763002 s. 272 ,Returning DM Positions 27.272540 s. 273 ,Returning DM Positions 36.897714 s. 274 ,Returning DM Positions 23.519835 s. 275 ,Returning DM Positions 30.785704 s. 276 ,Returning DM Positions 23.773792 s. 277 ,Returning DM Positions 31.145520 s. 278 ,Returning DM Positions 48.254341 s. 279 ,Returning DM Positions 316.176956 s. 280 ,Returning DM Positions 28.322149 s. 281 ,Returning DM Positions 33.514816 s. 282 ,Returning DM Positions 22.171488 s. 283 ,Returning DM Positions 73.418332 s. 284 ,Returning DM Positions 31.373158 s. 285 ,Returning DM Positions 29.498157 s. 286 ,Returning DM Positions 26.208930 s. 287 ,Returning DM Positions 34.666309 s. 288 ,Returning DM Positions 206.032375 s. 289 ,Returning DM Positions 23.002695 s. 290 ,Returning DM Positions 27.909436 s. 291 ,Returning DM Positions 22.943056 s. 292 ,Returning DM Positions 20.465738 s. 293 ,Returning DM Positions 22.692198 s. 294 ,Returning DM Positions 35.717975 s. 295 ,Returning DM Positions 33.996248 s. 296 ,Returning DM Positions 28.469246 s. 297 ,Returning DM Positions 23.901075 s. 298 ,Returning DM Positions 23.166638 s. 299 ,Returning DM Positions 22.924115 s. 300 ,Returning DM Positions 26.428412 s. 301 ,Returning DM Positions 21.613242 s. 302 ,Returning DM Positions 36.596315 s. 303 ,Returning DM Positions 23.771823 s. 304 ,Returning DM Positions 33.350181 s. 305 ,Returning DM Positions 27.001312 s. 306 ,Returning DM Positions 39.864418 s. 307 ,Returning DM Positions 22.848017 s. 308 ,Returning DM Positions 21.168089 s. 309 ,Returning DM Positions 40.674191 s. 310 ,Returning DM Positions 23.071757 s. 311 ,Returning DM Positions 23.319320 s. 312 ,Returning DM Positions 20.835019 s. 313 ,Returning DM Positions 36.229717 s. 314 ,Returning DM Positions 33.049161 s. 315 ,Returning DM Positions 26.786599 s. 316 ,Returning DM Positions 22.713673 s. 317 ,Returning DM Positions 30.654183 s. 318 ,Returning DM Positions 40.599154 s. 319 ,Returning DM Positions 36.636168 s. 320 ,Returning DM Positions 31.809886 s. 321 ,Returning DM Positions 33.322869 s. 322 ,Returning DM Positions 26.006976 s. 323 ,Returning DM Positions 23.468815 s. 324 ,Returning DM Positions 229.984027 s. 325 ,Returning DM Positions 41.038028 s. 326 ,Returning DM Positions 24.127290 s. 327 ,Returning DM Positions 145.922044 s. 328 ,Returning DM Positions 40.646327 s. 329 ,Returning DM Positions 47.194231 s. 330 ,Returning DM Positions 570.882853 s. 331 ,Returning DM Positions 21.435304 s. 332 ,Returning DM Positions 21.326541 s. 333 ,Returning DM Positions 26.414368 s. 334 ,Returning DM Positions 56.995868 s. 335 ,Returning DM Positions 28.150712 s. 336 ,Returning DM Positions 24.039218 s. 337 ,Returning DM Positions 38.550505 s. 338 ,Returning DM Positions 24.168785 s. 339 ,Returning DM Positions 21.817693 s. 340 ,Returning DM Positions 22.902155 s. 341 ,Returning DM Positions 25.154035 s. 342 ,Returning DM Positions 29.292105 s. 343 ,Returning DM Positions 28.569669 s. 344 ,Returning DM Positions 36.914976 s. 345 ,Returning DM Positions 22.288218 s. 346 ,Returning DM Positions 21.803264 s. 347 ,Returning DM Positions 21.774840 s. 348 ,Returning DM Positions 42.763051 s. 349 ,Returning DM Positions 22.517285 s. 350 ,Returning DM Positions 30.844080 s. 351 ,Returning DM Positions 21.505363 s. 352 ,Returning DM Positions 36.070710 s. 353 ,Returning DM Positions 23.194426 s. 354 ,Returning DM Positions 22.395653 s. 355 ,Returning DM Positions 27.906052 s. 356 ,Returning DM Positions 67.733488 s. 357 ,Returning DM Positions 22.617105 s. 358 ,Returning DM Positions 42.594370 s. 359 ,Returning DM Positions 29.905415 s. 360 ,Returning DM Positions 29.033982 s. 361 ,Returning DM Positions 34.716625 s. 362 ,Returning DM Positions 34.183171 s. 363 ,Returning DM Positions 30.334727 s. 364 ,Returning DM Positions 26.180832 s. 365 ,Returning DM Positions 26.093821 s. 366 ,Returning DM Positions 30.979216 s. 367 ,Returning DM Positions 29.474466 s. 368 ,Returning DM Positions 24.332565 s. 369 ,Returning DM Positions 22.197293 s. 370 ,Returning DM Positions 53.481345 s. 371 ,Returning DM Positions 178.519249 s. 372 ,Returning DM Positions 47.247952 s. 373 ,Returning DM Positions 25.966159 s. 374 ,Returning DM Positions 41.392763 s. 375 ,Returning DM Positions 39.934367 s. 376 ,Returning DM Positions 28.345356 s. 377 ,Returning DM Positions 46.177529 s. 378 ,Returning DM Positions 23.032810 s. 379 ,Returning DM Positions 34.932183 s. 380 ,Returning DM Positions 22.012015 s. 381 ,Returning DM Positions 113.736556 s. 382 ,Returning DM Positions 22.903738 s. 383 ,Returning DM Positions 25.592198 s. 384 ,Returning DM Positions 24.719515 s. 385 ,Returning DM Positions 21.909728 s. 386 ,Returning DM Positions 26.153578 s. 387 ,Returning DM Positions 23.789823 s. 388 ,Returning DM Positions 23.162080 s. 389 ,Returning DM Positions 26.753594 s. 390 ,Returning DM Positions 26.263075 s. 391 ,Returning DM Positions 101.164749 s. 392 ,Returning DM Positions 129.082483 s. 393 ,Returning DM Positions 41.936748 s. 394 ,Returning DM Positions 28.953769 s. 395 ,Returning DM Positions 22.590365 s. 396 ,Returning DM Positions 21.968985 s. 397 ,Returning DM Positions 28.682487 s. 398 ,Returning DM Positions 32.846702 s. 399 ,Returning DM Positions 22.329725 s. 400 ,Returning DM Positions 30.865709 s. 401 ,Returning DM Positions 28.868161 s. 402 ,Returning DM Positions 22.068793 s. 403 ,Returning DM Positions 25.523471 s. 404 ,Returning DM Positions 22.393050 s. 405 ,Returning DM Positions 23.983187 s. 406 ,Returning DM Positions 20.810259 s. 407 ,Returning DM Positions 25.810385 s. 408 ,Returning DM Positions 26.075842 s. 409 ,Returning DM Positions 25.674633 s. 410 ,Returning DM Positions 34.985300 s. 411 ,Returning DM Positions 29.858105 s. 412 ,Returning DM Positions 26.776467 s. 413 ,Returning DM Positions 21.324376 s. 414 ,Returning DM Positions 21.916645 s. 415 ,Returning DM Positions 36.489125 s. 416 ,Returning DM Positions 26.844995 s. 417 ,Returning DM Positions 23.138471 s. 418 ,Returning DM Positions 30.624277 s. 419 ,Returning DM Positions 28.346509 s. 420 ,Returning DM Positions 26.324284 s. 421 ,Returning DM Positions 24.490266 s. 422 ,Returning DM Positions 36.595900 s. 423 ,Returning DM Positions 35.991509 s. 424 ,Returning DM Positions 40.536067 s. 425 ,Returning DM Positions 41.070200 s. 426 ,Returning DM Positions 24.991530 s. 427 ,Returning DM Positions 40.494021 s. 428 ,Returning DM Positions 32.526759 s. 429 ,Returning DM Positions 24.764166 s. 430 ,Returning DM Positions 21.914904 s. 431 ,Returning DM Positions 24.232727 s. 432 ,Returning DM Positions 26.683221 s. 433 ,Returning DM Positions 23.850124 s. 434 ,Returning DM Positions 38.025230 s. 435 ,Returning DM Positions 23.363063 s. 436 ,Returning DM Positions 23.066884 s. 437 ,Returning DM Positions 436.106050 s. 438 ,Returning DM Positions 24.091832 s. 439 ,Returning DM Positions 25.066797 s. 440 ,Returning DM Positions 26.349757 s. 441 ,Returning DM Positions 20.335750 s. 442 ,Returning DM Positions 27.254616 s. 443 ,Returning DM Positions 23.144879 s. 444 ,Returning DM Positions 21.757001 s. 445 ,Returning DM Positions 43.739453 s. 446 ,Returning DM Positions 174.483205 s. 447 ,Returning DM Positions 27.688447 s. 448 ,Returning DM Positions 31.449643 s. 449 ,Returning DM Positions 27.063838 s. 450 ,Returning DM Positions 55.589213 s. 451 ,Returning DM Positions 56.283825 s. 452 ,Returning DM Positions 23.844368 s. 453 ,Returning DM Positions 23.604308 s. 454 ,Returning DM Positions 24.660988 s. 455 ,Returning DM Positions 65.463399 s. 456 ,Returning DM Positions 72.533486 s. 457 ,Returning DM Positions 23.385254 s. 458 ,Returning DM Positions 22.965261 s. 459 ,Returning DM Positions 22.739172 s. 460 ,Returning DM Positions 27.548026 s. 461 ,Returning DM Positions 25.431535 s. 462 ,Returning DM Positions 30.999227 s. 463 ,Returning DM Positions 29.952240 s. 464 ,Returning DM Positions 24.732376 s. 465 ,Returning DM Positions 25.449651 s. 466 ,Returning DM Positions 24.236293 s. 467 ,Returning DM Positions 103.359702 s. 468 ,Returning DM Positions 25.504343 s. 469 ,Returning DM Positions 23.619103 s. 470 ,Returning DM Positions 35.609036 s. 471 ,Returning DM Positions 30.573133 s. 472 ,Returning DM Positions 42.818583 s. 473 ,Returning DM Positions 24.456975 s. 474 ,Returning DM Positions 24.950212 s. 475 ,Returning DM Positions 26.776945 s. 476 ,Returning DM Positions 53.317369 s. 477 ,Returning DM Positions 36.714204 s. 478 ,Returning DM Positions 27.444448 s. 479 ,Returning DM Positions 22.642990 s. 480 ,Returning DM Positions 24.631790 s. 481 ,Returning DM Positions 35.226649 s. 482 ,Returning DM Positions 25.379335 s. 483 ,Returning DM Positions 32.527921 s. 484 ,Returning DM Positions 23.097038 s. 485 ,Returning DM Positions 113.051942 s. 486 ,Returning DM Positions 176.316620 s. 487 ,Returning DM Positions 24.510122 s. 488 ,Returning DM Positions 29.268195 s. 489 ,Returning DM Positions 24.624661 s. 490 ,Returning DM Positions 37.973792 s. 491 ,Returning DM Positions 31.971331 s. 492 ,Returning DM Positions 26.714006 s. 493 ,Returning DM Positions 34.880574 s. 494 ,Returning DM Positions 36.247390 s. 495 ,Returning DM Positions 22.645120 s. 496 ,Returning DM Positions 32.959379 s. 497 ,Returning DM Positions 34.934201 s. 498 ,Returning DM Positions 129.075725 s. 499 ,Returning DM Positions 22.007288 s. 500 ,Returning DM Positions 24.472511 s. 501 ,Returning DM Positions 23.716115 s. 502 ,Returning DM Positions 21.149309 s. 503 ,Returning DM Positions 35.102476 s. 504 ,Returning DM Positions 39.811308 s. 505 ,Returning DM Positions 34.504229 s. 506 ,Returning DM Positions 26.408525 s. 507 ,Returning DM Positions 36.067188 s. 508 ,Returning DM Positions 22.020084 s. 509 ,Returning DM Positions 36.392552 s. 510 ,Returning DM Positions 20.847379 s. 511 ,Returning DM Positions 23.411104 s.
In [159]:
hp.mollview(stackR, title="Mollview image RING")
hp.graticule()
0.0 180.0 -180.0 180.0 The interval between parallels is 30 deg -0.00'. The interval between meridians is 30 deg -0.00'.
In [166]:
hp.mollview(hp.sphtfunc.smoothing(stackR,0.3), title="Mollview image RING")
hp.graticule()
Sigma is 437.962991 arcmin (0.127398 rad) -> fwhm is 1031.324031 arcmin Sigma is 0.000000 arcmin (0.000000 rad) -> fwhm is 0.000000 arcmin 0.0 180.0 -180.0 180.0 The interval between parallels is 30 deg -0.00'. The interval between meridians is 30 deg -0.00'.
In [162]:
LMAX = 2*NSIDE
clR = hp.anafast(stackR, lmax=LMAX)
ellR = np.arange(len(clR))
In [169]:
plt.figure(figsize=(10, 5))
plt.plot(ellR, ellR * (ellR + 1) * clR)
#plt.semilogy(ellR, ellR * (ellR + 1) * clR)
#plt.plot(ell, ell * (ell + 1) * cl/2000)
plt.xlabel("$\ell$")
plt.ylabel("$\ell(\ell+1)C_{\ell}$")
plt.title('R = 4 x halfMassRadius')
plt.grid()
#plt.xlim(0,17)
In [176]:
plt.figure(figsize=(10, 5))
#plt.plot(ellR, ellR * (ellR + 1) * clR)
plt.semilogy(ellR, clR)
#plt.plot(ell, ell * (ell + 1) * cl/2000)
plt.xlabel("$\ell$")
plt.ylabel("$C_{\ell}$")
plt.title('R = 4 x halfMassRadius')
plt.grid()
#plt.xlim(0,17)
END¶
In [177]:
start = time.time()
#velocities = readsnap(infile, 'vel', 'dm')
end = time.time()
print("Snapshot readed in %f s." %(end-start))
Returning DM Positions Snapshot readed in 0.263317 s.
In [188]:
from mpl_toolkits.mplot3d import * #3D-s ábrák alcsomagja
from ipywidgets import * #interaktivitáshoz szükséges függvények
In [195]:
plt.plot(stack,'.')
Out[195]:
[<matplotlib.lines.Line2D at 0x7f82adb576a0>]
In [255]:
ss = hp.sphtfunc.smoothing(stack,0.4)
xyzDist = hp.pix2vec(NSIDE, range(len(ss))) * (ss-0.95*ss.mean())
Sigma is 583.950655 arcmin (0.169864 rad) -> fwhm is 1375.098708 arcmin Sigma is 0.000000 arcmin (0.000000 rad) -> fwhm is 0.000000 arcmin
In [256]:
from plotly import *
from plotly.offline import *
from plotly.graph_objs import *
init_notebook_mode()
trace = Scatter3d(x=xyzDist[0,:], y=xyzDist[1,:], z=xyzDist[2,:], mode='markers',
marker=dict(
size=1,
colorscale='Viridis', # choose a colorscale
opacity=0.8
)
) # Adatok csoportosÃtása trace-be
adatok = [trace]
layout3d=Layout(title='3d')
fig_3D = Figure(data=adatok,layout=layout3d) # ábra objektum létrehozása
iplot(fig_3D) # ábra megjelenÃtése
In [262]:
plt.figure(figsize=(10,10))
plt.plot(xyzDist[0,:],xyzDist[1,:],'.')
Out[262]:
[<matplotlib.lines.Line2D at 0x7f82ade469b0>]
In [ ]: